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Abstract 

We study a versatile structurally favorable periodic -bonded carbon atomic planar sheet with C/^^ sym- 
metry by means of the first-principles calculations. This carbon allotrope is composed of carbon octagons 
and squares with two bond lengths and is thus dubbed as octagraphene. It is a semimetal with the Fermi 
surface consisting of one hole and one electron pocket, whose low-energy physics can be well described 
by a tight-binding model of vr-electrons. Its Young's modulus, breaking strength and Poisson's ratio are 
obtained to be 306 N/m, 34 A N/m and 0.13, respectively, which are close to those of graphene. The 
novel sawtooth and armchair carbon nanotubes as well as unconventional fullerenes can also be constructed 
from octagraphene. It is found that the Ti-absorbed octagraphene can be allowed for hydrogen storage with 
capacity around 7.76 wt%. 
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I. INTRODUCTION 



n 

The discovery of graphene [1] opens a new era of fundamental research and technological 
applications of carbon-based materials. Graphene is an elegant two-dimensional (2D) material 
exhibiting exceptional electrical, thermal, mechanical and optical properties, and its low-energy 
behavior of electrons is described by a (2+l)-dimensional relativistic quantum theory tA^- Stim- 
ulated by the tremendous interest of graphene, a great number of novel carbon structures were 



proposed and investigated in the past years (e.g. Refs. p MlOh ). Among others, two intriguing 
sp-sp'^ -hondtd monolayers of elemental carbon, graphyne [U] and graphdiyne Q, [isl], and an 



insulating 2D hydrocarbon layer, graphane [150, are typical examples, all of which were ex- 
perimentally obtained. It appears that there is a growingly extensive interest to seek for more 2D 
allotropes of carbon with amazing physical and chemical properties as well as possibly diverse 
applications. 

Graphene has a perfect honeycomb lattice with C^y symmetry. Recent studies [|6|, 



170 show 



that the line defects with octagons and pentagons could be formed self-assembly in graphene. The 
formation of small carbon molecules with various polygon rings has also been studied from exper- 
imental aspects [18]. Inspired by these observations, we propose a stable 2D periodic atomic sheet 
consisting of carbon octagons, coined as octagraphene, which possesses intriguing properties and 
might be synthesized experimentally by line defects or acetylene scaffolding ways. As shown 
in Fig. 1(a), octagraphene comprises of carbon octagons and squares with two bond lengths that 
forms a square lattice with C^^ symmetry. By means of the first-principles calculations, it is shown 
that octagraphene is energetically and kinetically stable, and is more favorable in energy than gra- 
phyne and graphdiyne, although it is metastable against graphene. By rolling octagraphene along 
specific directions, we can get novel sawtooth and armchair carbon nanotubes [Fig. 1(b)] that 
are metallic, which can escape from the difficult separation of usual metallic and semiconducting 
carbon nanotubes uM \ by wrapping octagraphene and making use of pentagons, hexagons or hep- 
tagons as caps, many stable unconventional fuUerenes such as C36, C48, Cso, Cge and so forth [Fig. 
1(c)] can be constructed. The electronic properties of octagraphene manifest that it is a semimetal 
with small density of states (DOS) around the Fermi level, whose Fermi surface consists of one 
hole and one electron pockets centered at V and M points, respectively. The low-energy physics 
of octagraphene can be well described by a tight-binding model of tt electrons. Octagraphene has 
the density of 0.68 mg/m?. Young's modulus of 306 N/m, breaking strength of 34.4 N/m and 



Poisson's ratio of 0.13. It could be the strongest periodic carbon sheet after graphene till now. 
An energy gap can be opened by doping boron nitrogen pairs. The possible routes for obtaining 
octagraphene are suggested. The Ti-absorbed octagraphene can also be used for hydrogen storage 
with capacity of 7.76 wt%. With these amazing characters, we anticipate that octagraphene, sim- 
ilar to its cousin — graphene, once being obtained, will likewise have a crucial impact in physics, 
chemistry, materials and information sciences. 

The rest of this article is organized as follows. In Sec. II, the calculational methods are de- 
scribed. In Sec. Ill, the geometrical structure, mechanical properties, electronic structures and 
low-energy physics of octagraphene are presented. In Sec. IV, new carbon nanotubes and uncon- 
ventional fuUerenes from octagraphene are given. The band engineering and hydrogen storage of 
octagraphene are discussed in Sec. V. Finally, a conclusion is briefly presented. 



II. CALCULATIONAL METHODS 



Both the first-principles calculations within the framework of the density-functional theory 



(DFT) 



2 If] and the tight-binding approximation (TBA) [2^ are invoked to study this 2D sys- 



tem. For the first-principles method, primary calculations were performed within the Vienna ab 
initio simulation package (VASP) [|23|, with the projector augmented wave (PAW) method 
ll25n . Both the local density approximation (LDA) in the form of Perdew-Zunger ll26n and gener- 
alized gradient approximation (GGA) developed by Perdew and Wang ll27n were adopted for the 
exchange correlation potentials. The plane- wave cutoff energy is taken as 400 eV. The Monkhorst- 
Pack scheme for k-point samplings with 32x32x 1 mesh for planar sheets and 1 x 1 x 100 mesh for 
nanotubes were used to sample the Brillouin zone [28]. The supercells are used for calculations of 
isolated sheet structures, and the distance between two layers is about 10 A to avoid interactions. 
The geometries were optimized when the remanent Hellmann-Feynman forces on the ions are less 
than 0.01 eV/A. Phonon calculations with a k-mesh of 8x8x1 are performed using the density 
functional perturbation theory [29.1 . 



(a) octagraphene (b) OCNT(10,0) OCNT(8,8) 




FIG. 1. (Color online) Schematic depiction of the structures of (a) octagraphene, where a unit cell is indi- 
cated with the unit vectors ai and a2; (b) single-wall sawtooth (10,0) and armchair (8,8) carbon nanotubes 
rolled from octagraphene; and (c) unconventional fuUerenes of C36, C48, Cso and Cqq obtained from octa- 
graphene. 

III. OCTAGRAPHENE 

A. Geometrical structure and mechanical properties 

The optimized geometric structure of octagraphene is presented in Fig. 1(a). Carbon atoms 
in octagraphene form a square lattice with the lattice constant qq = 3.45 A, where a unit cell 
contains four carbon atoms. Similar to graphene, every carbon atom in octagraphene has three 
sp^ -bonded nearest neighbors, forming three a bonds. In addition, octagraphene has two types of 
bond lengths, 1.48 A of intra-squares and 1.35 A of inter- squares, and two bond angles 90° and 
135°. 

The cohesive energy Ec per atom as a function of lattice constant (a/ao) is given in Fig. 2(a) 
for graphene, octagraphene, graphyne and graphdiyne, respectively. It is clear that a single mini- 
mum of —Ec appears at a/ao = 1 for these 2D structures, suggesting that they are energetically 
stable. Note that the minimum of —Ec for octagraphene is smaller than those of graphyne and 
graphdiyne but larger than that of graphene, suggesting octagraphene is more energetically sta- 
ble than graphyne and graphdiyne but metastable against graphene. Fig. 2(b) gives the phonon 
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FIG. 2. (Color online) (a) Cohesive energy per atom as a function of lattice constant a for graphene, 
octagraphene, graphyne and graphdiyne, respectively, where ao is the optimized lattice constant, (b) Phonon 
spectra of octagraphene. 

spectra of octagraphene. No imaginary phonon modes are found, implying that octagraphene is 
also kinetically stable. We performed quantum molecular dynamics simulations and found that at 
temperature 500 K the planar structure of octagraphene still maintains. Considering that graphyne 
and graphdiyne wtrt already obtained experimentally, it is reasonable to believe that octagraphene 
can be synthesized in lab or likely already exists in Nature. In fact, carbon rings w^ith four or eight 



sides have already been observed in some experiments ml 

For a comparison, the cohesive energy Ec, density p, bond length Icc^ energy gap Eg, Young's 
modulus E, breaking strength a and Poisson's ratio u for graphene, octagraphene, graphyne and 
graphdiyne, respectively, are collected in Table I. Notice that apart from the experimental data of 
graphene taken from literature, all data presented here wtre calculated by ourselves using the same 
method. As octagons in octagraphene are larger in size than hexagons in graphene, octagraphene 
has the density (0.68 mg/rri^) slightly smaller than that of graphene (0.77 mg/rri^). 

Nov^ let us discuss the mechanical properties of octagraphene. As it has a C/^^ symmetry, only 
three elastic constants cij are independent. Our DFT calculations estimate that the values of cn, 
Ci2 and C44 of octagraphene are 296, 46, and 49 N/m, respectively. The Young's modulus E 
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FIG. 3. (Color online) (a) Stress-strain curves of graphene (GE), octagraphene (OGE), graphyne (GY) and 
graphdiyne (GDY) along armchair (AC) and zigzag (ZZ) or sawtooth (ST) directions. The experimental 



data of graphene are taken from OOll . (b) Valence electron density of octagraphene. 



of octagraphene is obtained to be 306 and 168 N/m along the sawtooth and armchair directions 
(see below), respectively, in comparison to those of graphene 350 N/m, graphyne 245 N/m 
and graphdiyne 174 N/m. The stress-strain curves for graphene, octagraphene, graphyne and 
graphdiyne under a uniaxial tension along two non-chiral directions are illustrated in Fig. 3(a). 
The ultimate tensile strength or the breaking strength a of octagraphene is 34.4 and 27.3 N/m 
along the sawtooth and armchair directions, respectively, which are quite comparable with those 
of graphene 38.4 A^/m (zigzag direction) and 34.8 N/m (armchair direction) but larger than those 
of graphyne (23.4 and 19.5 N/m) and graphdiyne (16.8, 13.4 N/m). The experimental values 
OOn of E and a for graphene are 340 ± 50 N/m and 42 ± 4 N/m, respectively, showing that 
our DFT calculations are quite reliable. The Poisson's ratio of octagraphene is 0.13 (sawtooth 
direction) and 0.47 (armchair direction), which are comparable with 0.21 (zigzag direction) and 
0.17 (armchair direction) of graphene. The Poisson's ratio of graphyne and graphdiyne can be 
found in Table I. From above results we see that octagraphene has the mechanical properties very 
similar to graphene, and might be the strongest carbon atomic sheet after graphene till now. 



B. Electronic structures 



The valence electron density, electronic structures and the DOS of octagraphene are obtained 
by means of the DFT and TBA calculations, respectively. From Fig. 3(b) one may see that more 
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TABLE I. The symmetry, lattice constant (ao), bond lengths Icc^ plane density p, cohesive energy (Ec), 
energy gap Eg between the bottom of conduction band and the top of valence band, Young's modulus E, 
breaking strength a and Poission's ratio u for graphene, octagraphene, graphyne and graphdiyne, respec- 
tively. 



Structure method symmetry ao (A) Icc (A) p {mg/m^) Ec (eV/atom) Eg (eV) E (N/m) a (N/m) v 
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FIG. 4. (Color online) (a) Electronic band structures and density of states (DOS) of octagraphene calculated 
by the density functional theory (DET) and the tight-binding approximation (TBA) methods, (b) Eermi 
surface of octagraphene, where the hole packet (red) is at F point, and the electron pocket (blue) is at M 
point, (c) Energy bands of octagraphene nanotubes (3,3), (5,5), (12,12), (5,0), (8,0) and (16,0). 
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electrons are distributed on the inter-square bonds than on the intra-square bonds in octagraphene, 
which is well consistent with its structural feature, showing an unequal sp^ hybridizing character. 
The electronic structures in Fig. 4(a) show that octagraphene is a semimetal, as its top of valence 
band and the bottom of conduction band are located at V and M points in the Brillouin zone, 
respectively, giving rise to an indirect negative band gap (-4.03 eV), and leading to the Fermi level 
passes through both the conduction and valence bands. Thus, in contrast to graphene, the charge 
carriers in octagraphene have both holes and electrons around the Fermi level. Note that graphene 
is also a semimetal or a zero-gap semiconductor because its valence and conduction bands touch 
at K points, while graphyne and graphdiyne are semiconductors with narrow band gap of 0.48 and 
0.51 eV, respectively. 

The projected DOS of octagraphene shows that the energy bands near the Fermi surface are 
predominantly contributed by the 2pz orbital of carbon atom, forming the same 7r-bond as that 
in graphene. The Fermi surface of octagraphene is plotted in Fig. 4(b), which consists of a hole 
pocket at V point and an electron pocket at M point. Since each unit cell contributes four electrons 
with two bonding states and two antibonding states, each point in the Brillouin zone should be 
filled with 2x2 electrons, and the number of unoccupied states in valence band is just equal to 
that of occupied states in conduction band, resulting in that the hole and electron pockets share the 
same area. 



C. Low-energy physics 

To describe the low-energy physics near the Fermi surface of octagraphene, we propose a tight- 
binding model which contains only tt electrons with nearest neighbor hoppings. The Hamiltonian 
can be written as 

where and cj^ are the electron annihilation and creation operators at site i with spin cr, (i, j) 
runs over the nearest neighbors within the same square and t is the corresponding hopping ampli- 
tude of electrons, <^ i,j ^ runs over the nearest neighbors connecting different squares and f is 
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the corresponding hopping amphtude. In the matrix form, the Hamiltonian reads 

/ t fe-'^y t \ 
t t t'e'^- 



t'e'^y t t 



(2) 



\ t t'e-'^^ t J 

where kx and ky are x and y components of momentum of tt electrons, and the lattice spacing 
ao = 1 is taken. The dispersion E(k) of octagraphene within the TBA is determined by 



E{kY - 2{t'^ + 2t^)E{kf + 4t't2(cos k^ + cos ky)E{k) 
+ - (cos kx + cos ky)] = 0. 



(3) 



The TBA results are in good agreement with those obtained by the DFT calculations [Fig. 4(a)], 
implying that the low-energy properties of tt electrons in octagraphene can be well described by 
the TBA Hamiltonian, where the fittings give t = —2.5 eV and f = —2.9 eV. For graphene, 
t = -2.7 eV [341. 



D. Fabricating suggestions 

To obtain the structure of octagraphene experimentally, we suggest that one may make the line 
defects periodically in graphene, and then connect properly the broken bonds of carbon atoms, 
as demonstrated in Fig. 5. In addition, it may be possible to synthesize it through the acetylene 



scaffolding ways 



m 



IV. OCTAGRAPHENE NANOTUBES AND FULLERENES 

By rolling octagraphene, the novel single-walled carbon nanotubes can also be obtained. For 
this purpose, let us define the chiral vector of octagraphene carbon nanotube (OCNT) by Ch = 
nai + ma2 = {n,m), where n and m are integers with < \m\ < n because of the C^y symmetry. 
We call (n, 0) the sawtooth OCNT (m = 0, 6> = 0°), and (n, n) the armchair OCNT (m = n, 
9 = 45°), where the chiral angle is defined as cos 9 = nj\/n^ + m^. The tube diameter is given by 

= IC^^I/tt = a^\/n^ + m? /ti. The geometric structures of sawtooth OCNT (10,0) and armchair 
OCNT (8,8) are presented in Fig. 1(b) as examples, whose structural stabilities were verified by 
the DFT calculations. The energy bands of OCNT (3,3), (5,5), (12,12), (5,0), (8,0) and (16,0) are 
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FIG. 5. (Color online) A schematic demonstration how to utilize the line defects in graphene to obtain 
octagraphene. 



given in Fig. 4(c), which show that regardless of the curvature effect, the single-walled OCNTs, 
either sawtooth or armchair, are all metallic. 

In addition, by wrapping octagraphene and using pentagons, hexagons or heptagons as caps, 
many structurally stable unconventional carbon fuUerenes can be built. The optimized examples 
like C36, C48, Cso and Cge are shown in Fig. 1(c). These results show that octagraphene is indeed 
a versatile 2D crystalline allotrope of carbon after graphene. 

V. BAND ENGINEERING AND HYDROGEN STORAGE 
A. Boron nitrogen pair doping 



Similar to graphene and carbon nanotubes[135Ll36[], the semimetallic octagraphene may be ap- 
plicable in the band engineering via the chemical doping. In fact our calculation reveals that a 
band gap opens if the octagraphene is substitutional^ doped with boron and nitrogen (B-N) pairs, 
as shown in Fig. [6l where eleven structures for different doping configurations and a pure B-N 
octa-square structure are considered. All these doped structures are energetically favorable. In Fig. 
[6] (a), the structures of number 1-11 are those of octagraphene of 4x4 supercell with substituting 
carbon (blue ball) by boron (red ball) and nitrogen (green ball) atom pairs, and the structure of 
number 12 is the pure B-N octagon- square structure of 3x3 supercell. In the first structure two 
carbon neighbors are replaced by one B-N pair; in the second (third) structure a square (octagon) 
is substituted by two (four) B-N pairs; in the fourth (sixth, seventh or ninth) one a nanoribbon in 
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FIG. 6. (a) Schematic structures of boron nitrogen pairs substitutionally doped octagraphene with different 
doping configurations, where carbon, boron, and nitrogen atoms are colored in blue, red, and green, respec- 
tively, (b) The corresponding band gaps (E^) for different structures of (a), where the doping concentration 
is from small to large in sequence. 
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FIG. 7. (Color online) There are five possible adsorption sites of metal atoms on octagraphene, including 
one top site on C atoms (T), two bridge sites on top of intra-square (B) and inter-square (B^) C-C bond, 
respectively, and two hollow sites on top of square (H) and octagon (H'), respectively. 



octagraphene is substituted by four (eight, twelve or sixteen) B-N pairs; in the fifth one a 2x2 
supercell is replaced by eight B-N pairs; in the eighth one a cross is replaced by fourteen B-N 
pairs; the tenth one is the dual structure of the eighth; the eleventh is the dual structure of the fifth. 
The band gap varies significantly for different doping configurations, as shown in Fig. [6](b). It is 
seen that the boron nitrogen pair doped octagraphene can open a gap from 0.009 eV (structure 2) 
to 1.866 eV (structure 11), which can gain some insights for nanoelectronics. It is noting that the 
band gaps of most of structures considered here are around 0.2 eV, even when the doping concen- 
tration is as large as 50.82% (structure 9), which is quite different from the doped graphene [1371] . 
Considering that the GGA usually underestimates the band gap, the present doping study gives 
useful information for the band engineering in octagraphene. 



B. Hydrogen storage 

Next, let us consider the possibility of hydrogen storage in octagraphene. In general, there are 
two ways to realize this purpose. One is that hydrogen atoms can be absorbed on carbons via the 



chemical bond to form a new hydrocarbon sheet, similar to graphane u4\ : the other is that the 
hydrogen molecules can be adsorbed physically on some metal atoms (e.g. Ti, Ca, Al, etc.) that 
are deposited on the sheet of octagraphene. In the former case, the hydrogen storage capacity of 
octagraphene can reach 7.7 wt% as graphane does In the latter case, the hydrogen capac- 
ity depends on the adsorption configuration of metal atoms. There are five possible adsorption 
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FIG. 8. (Color online) The configuration of four hydrogen molecules adsorbed on each Ti atom at Hi 
site, (a)-(d) reveal the single-side hydrogen storage configurations with one, two, three and four hydrogen 
molecules adsorbed on each Ti atom, respectively; (e)-(h) show the double-side hydrogen storage configu- 
rations with one, two, three and four H2 adsorbed on each Ti atom, respectively; where C, Ti and H atoms 
are colored in dark blue, bluish green and brick red, respectively. 

positions for metal atoms (see Fig. [7]), including one top site on C atoms (T), two bridge sites 
on top of intra-square (B) and inter-square (B^) C-C bond, respectively, and two hollow sites on 
top of square (H) and octagon (H^), respectively. To analyze the stability of these five adsorption 
configurations, we calculate the cohesive energy E'^-m with M the adsorbed metal atom by 



where Eocta^ Em and EM-octa are the total energies of pure octagraphene, M atom and octa- 
graphene with adsorbed metal M, respectively. For Ti doping (M=Ti), H site is the most stable 
position with cohesive energy of 2.71 eV per Ti atom. 

Figs. [8ta)-(d) present single-side hydrogen storage configurations with one, two, three and four 
hydrogen molecules adsorbed on each Ti atom. The average adsorption energy per H2 is 0.207 
eV, 0.361 eV, 0.336 eV and 0.325 eV, respectively. Figs. [8] (e)-(h) give double-side hydrogen 
storage configurations with one, two, three and four H2 adsorbed on each Ti atom. The average 
adsorption energy per H2 is 0.235 eV, 0.355 eV, 0.336 eV and 0.325 eV, respectively. In this case, 
the hydrogen storage capacity could reach 7.76 wt% for Ti [as shown in Fig. [8][h)], which is very 



Eb-M — [Eocta + — EM-octa] /"^^ 



(4) 
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close to that of Ca- adsorbed graphene, say 8.4 wt% ll4Qh . We also perform the same calculation for 
Al and Ca atoms (M=Al or Ca), respectively, and found that the adsorption energy of hydrogen 
molecule is too small. The adsorption energy of Ti-adsorbed octagraphene is quite promising, say, 
about 0.3 eV/H2. It appears that the hydrogen storage capacity of octagraphene is close to that of 
fuUerenes [38], carbon nanotubes [39] and graphene I40I]. 



VI. CONCLUSION 



In conclusion, we propose a periodic carbon sheet by the first-principles calculations. This 
novel carbon monolayer is a -hybridized structure and coined as octagraphene. Unlike 
graphene, octagraphene is composed of carbon squares and octagons possessing C^y symme- 
try, and is more favorable in energy than graphyne and graphdiyne. Its optimized structure has 
two bond lengths with 1.48 A within each square and 1.35 A between adjacent squares. The band 
structure shows that it is a semimetal with Fermi surface consisting of one hole pocket and one 
electron pocket. In order to describe the low-energy physics of octagraphene, we also suggest a 
tight binding model of tt electrons which is uncovered to agree well with DFT calculations, and is 
quite helpful for further exploring the properties of octagraphene in a magnetic field or emergent 
quantum phenomena such as quantum Hall effect, Kondo effect, and so on in octagraphene. The 
intriguing mechanical properties of octagraphene are also obtained, which are observed close to 
those of graphene. Similar to graphene, we find that octagraphene can also be rolled seamlessly 
to form stable carbon nanotubes and unconventional fuUerenes. By substitutionally doping boron 
nitrogen pairs, we disclose that the semimetallic octagraphene can be turned into a semiconductor 
with a band gap, the property will be useful for possible applications in nanoelectronics. This new 
carbon allotrope has also the possibility for hydrogen storage. Our calculation shows that for the 
2D octagraphene sheet the hydrogen storage capacity could reach 7.76% via Ti-adsorption that is 
close to that of graphene. Finally, we would like to mention that one might obtain this new kind 
of carbon atomic sheet via making line defects in graphene or through the acetylene scaffolding 
ways. 
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